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APPROXIMATE TREATMENT OF V/STOL WALL INTERFERENCE 
FOR CLOSED CIRCULAR TUNNELS 

By Harry H. Heyson 
Langley Research Center 

SUMMARY 

An approximate treatment of V/STOL wall interference in a circular tunnel indi- 
cates that the interference factors at the model for this tunnel should be of the same 
order of magnitude as those presently available for a square tunnel of equal cross- 
sectional area. There is a greater degree of uncertainty with respect to the lateral and 
longitudinal distributions of interference; however, the available results for the square 
tunnel should be reasonably close to those of the circular tunnel provided that the model 
is relatively small in comparison with the test section. 

INTRODUCTION 

References 1 to 5 have provided theoretical treatments of wind-tunnel interference 
for V/STOL models in rectangular tunnels. As yet, however, there has been no equiva- 
lent treatment of such models in tunnels of circular cross section despite the fact that 
several such turmels are regularly used for V/STOL studies both in this country and 
abroad. 

References 1 to 5 obtained interference factors by the use of an external image 
system to represent the flat boundaries of the tunnels. In the case of a circular tunnel, 
only limited information can be obtained in this manner (for example, ref. 6); however, 
other techniques (refs. 7 and 8) can be used for solution of the complete three-dimensional 
problem. A start in this direction was made recently by Michel (ref. 9) who obtained the 
interference for a single doublet of arbitrary inclination and location within a circular 
tunnel; however, the study of reference 9 was never completed as the integrations 
required to extend the results to a complete wake were not carried out. 

The present study does not attempt a rigorous treatment of the circular tunnel. 
Instead, two square sections of different orientation, the interference in which should 
bound the interference in the circular tunnel, are examined. The computed results 
should indicate the magnitude of the difference to be expected between circular and square 
test sections. In view of the approximate nature of the analysis, only centered vanishingly 
small models are considered herein. 


SYMBOLS 


Because of the limitations of the computer-controlled plotter used to prepare the 
figures, small variations in these symbols appear in the figures. 

Ajn momentum area of lifting system 

Ax test-section cross-sectional area 



Z X 

Ao = Ai + — cos Y sin v 

^ ^ D D 

Cl lift coefficient, Lift/qS 

D semidiagonal length of square cross section 

H semiheight of square tunnel 

K function related to induced velocities of wake in free air 

m,n integers 

q dynamic pressure 

R radius of tunnel with circular cross section 

S wing area 

t integer, equal to zero for lift forces and 1 for drag forces 

u,v,w induced velocitie.s directed parallel to X-, Y-, and Z-axis, respectively 

Uq mean value of longitudinal Induced velocity 

wq mean value of vertical induced velocity 

x,y,z distances measured along the X-, Y-, and Z-axis, respectively 
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X,Y,Z 


6 


Aa 

X 


Cartesian coordinate axes centered in model and tunnel; X-axis positive 
rearward, Z-axis positive upward, and the Y-axis oriented to form a 
right-hand coordinate system 


wind-tunnel interference factor defined in terms of momentum area and 

^m j 

mean induced velocity; for example, Auj_, = l ~ z — ^0 

■^Itl 

= 6w,d uq 


change in angle of attack caused by wall interference, radians 

wake skew angle, angle at which wake leaves the model, measured positive 
rearward from negative Z-axis to center of wake 


Subscripts: 

u pertaining to longitudinal velocities 

V pertaining to lateral velocities 

w pertaining to vertical velocities 

D pertaining to drag forces 

L pertaining to lift forces 

°o pertaining to a semi -infinite wake 

THEORY 
General Approach 

The initial solution for a wind-tunnel interference problem was for the circular 
tunnel with a model having an undeflected wake (ref. 6). After obtaining the interference 
for this case, Prandtl noted that the interference in a square tunnel should be about the 
same as that in a circular tunnel of equal cross-sectional area. This assumption was 
later borne out when Glauert (ref. 10) obtained an interference factor /defined in this case 
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AqA't\ 

as 6 = of 0.137 for a small model in the square tunnel as compared with the cor- 

SCl/ 

responding value of 0.125 obtained by Prandtl in the circular tunnel. 

The situation becomes more complicated when the wake is deflected substantially 
downward from the horizontal. References 1 to 5 have indicated that for large deflections 
the interference factors increase and, furthermore, that the increase in interference is 
largely due to a greatly increased effect of the floor upon the phenomenon. Thus, it might 
be expected that the "floor" of the circular tunnel, being concave upward, would further 
increase the interference over that for the square tunnel. On the other hand, if the tun- 
nels are of equal cross-sectional area (fig. 1), the lower regions of the circular tunnel 
are somewhat farther below the model than the floor of the equivalent square tunnel. The 
result of this increased distance should tend to counter to some degree the concavity of 
the circular tunnel. 

If the square tunnel is now rotated 45° (to form what will be termed herein a dia- 
mond tunnel), it will be seen that the differences between the two square configurations 
are greater than the differences between the square and circular configurations. 

(Observe that the diamond section has a greater effective concavity of the floor than the 
circular configuration and that the center of the floor is farther from the model in the 
diamond section.) The advantage of considering the problem in this light is that the 
interference in both the square and diamond sections can be calculated from image sys- 
tems with little mathematical difficulty. Once the interference factors for both rectangu- 
lar configurations are available, it is reasonable to assume that the interference in the 
circular tunnel should lie within these factors. For convenience, the present study takes 
the interference factors for the circular tunnel to be the arithmetic mean of those for the 
square and diamond sections. 


The Deflected Wake 

The wake assumed herein is identical to that of reference 3. It originates at the 
model, passes downward and rearward in a straight line defined by the angle x mea- 
sured from the vertical tunnel axis. In free air, the wake continues to infinity. (See 
fig. 2(a).) In the tunnel, however, it meets the lower boundary. At this point (as in the 
upper half of fig. 2(b)), the wake is assumed to turn and flow off along the lower boundary. 
Although no real wake could behave in such a manner, the wake comprised of the two 
linear paths described appears to be a reasonable approximation to the actual curved path 
of the wake, at least for conditions which do not result in completely gross distortions of 
flow throughout the entire tunnel. (See refs. 11 and 12.) 

The wake of the vanishingly small model considered herein is a simple doublet line 
wherein the axes of the doublets are inclined according to the induced drag-lift ratio of 
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the model. For convenience, the calculations are performed twice: once with the 
doublet axes vertical corresponding to pure lift forces; and once with the doublet axes 
horizontal corresponding to pure drag forces. Any arbitrary drag-lift ratio (or doublet 
inclination) can then be obtained by an appropriate superposition of the two results. This 
model is equally applicable to wings or more complicated devices such as rotors. (See 
ref. 13.) 


Square Section 

For the square test section, the analysis is identical to that of reference 3 and the 
interference factors can be obtained directly from that paper. As a computational con- 
venience, however, the values given herein were obtained by the use of a subroutine in 
the computer program given in the appendix. This subroutine (square) is merely a spe- 
cialization of the program given as appendix A of reference 5. 


Diamond Section 

Induced field in free air .- As noted previously, the wake assumed for the diamond 
tunnel is the same as that of reference 3. In the present usage, it is more convenient to 
nondimensionalize the expressions for the free-air wake in terms of the semidiagonal 
distance D across the tunnel. Thus for the centered model, the vertical induced veloc- 
ity due to lift is 



and, as a special case, when y = 90°, 


Kw 


2L L. £ 

D’D’D 


X = 90° 




(Ic) 


Similarly, the longitudinal induced velocity due to lift becomes 


Uoo = Wq 


-^m 

A>p 


X ^ ix y z 
■ D’D’D 


(2a) 
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where 




55 (fe * ^1 °°° ^ X 


Ai^A2 


AiW 


(2b) 


and, as a special case, when x = 90° 


K, 


■u 


lx = 90° 


^y_ A ^ 

D’D’Dj-Aj3 


(2c) 


In the present derivation, the lateral induced velocities are also required. Thus, 
differentiating the free-air potential of the semi-infinite wake of vertical doublets with 
respect to y and then nondimensionaUzing with respect to D 3delds 


Voo = Wf 


Aq* 


Ik (ILLL 

TT nO’D’D 


(3a) 


where 


K P y z\ _ 

Md’d’dJ ■ 


D! 


D 


D 


+ Aj cos X 1 






(3b) 


and, as a special case, when x = 90°, 




X ^ z 
D’D’b 


D D 


,2Ai - ^ 
3. 2\ ^ D 


X = 90° 


Ai A2 


(3c) 


The corresponding results for the wake of horizontal doublets are: The vertical 
induced velocity due to drag is 


Woo = UQ 




1^ /x y z^ 


(4a) 


where 


K, 


w 


/x ^ 
^D’D’D/ 


%% (§ + Ai cos x)(§ - Ai sin x 


Ai^A2 


Ai2a2^ 


(4b) 
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and, as a special case, when x = 90°, 


K 


w 


lD’D’d) 


= 90° 


D 

r 

Af 


(4c) 


The longitudinal induced velocity due to drag is 


u«o = Ur 


A»p 


Ik /2L L ^ 

^ ^D’D’D 


where 


K ( z 

Md’d’d 


ii ) ^{ i ) 


Ai^A2 


A1A2 


and, as a special case, when x = 9^^? 

X 


K, 


u 


^ Z. z= D 
D’D’Dr^^c 


(5a) 


(5b) 


(5c) 


lx = 90° 

The lateral induced velocity is obtained in a manner similar to equations (3) to 


obtain 


/X y z\ „ -^m 


Ik y. ^ 

■ 17 


(6a) 


where 


kI^L^- y I p ■ d 
A^’P’d] - ‘pIaiSaz 


- ^ - Ai sin x\ 


A 2. 2 

Ai A2 


(6b) 


and, as a special case, when x = 90°, 


Kx 


X y z\ _ D 


X = 90° 


D’d’D/ 


(6c) 
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Wake in the diamond test section .- Equations (1) to (6) proAdde only the induced 
velocity field of the semi -infinite wake in free air (fig. 2(a)). The field of the system in 
simple ground effect, under the present assumptions, can be obtained from the wake and 
image system indicated in figure 2(b). The field of the finite segment of inclined wake 
may be obtained by summing the fields of two semi-infinite wakes; one starting from the 
origin; the other, of equal but opposite strength, starting from the intersection of the 
original wake and the ground. The trailing leg of the wake also originates at the inter- 
section with the ground and always has a skew angle of 90°. 

The images below the ground, which are required in order to satisfy the condition 
of zero normal flow at the boundary, can be obtained in the same manner with suitable 
transpositions of the origin. Observe that a downwash with respect to the coordinate 
system of the image becomes an upwash with respect to the coordinate system of the real 
wake. 

The wake and image pattern required to produce zero normal flow through the two 
planes corresponding to the bottom corner of the diamond tunnel is somewhat more com- 
plicated. (See fig. 2(c).) Observe that the symmetries required to meet the boundary 
conditions in the corner necessitate three images. Of these three images, the first is 
identical to the image in the ground-effect system. The remaining two images lie in a 
horizontal plane through the corner. Because of the 90°-rotation of these images with 
respect to the real wake, it is the lateral induced velocities of the images which contribute 
to the vertical induced velocity field of the real wake. 

This fundamental cell of four wakes must be repeated in a doubly infinite pattern 
(as in fig. 2(d)) to satisfy conditions of zero normal flow through all four walls of the 
diamond section. Note that the basic unit of four cells repeats at intervals of 2(m - n)D 
in the Y-direction and at intervals of 2(m + n)D in the Z-direction. Thus, by defining 
the interference factors as in reference 3, the required superpositions yield 




6w.- = -:;l 


y - 2(m - n),| - 2<m + n)| + )> ^’n ‘ ‘ ‘ 2(m + n) + ll 

m=-«? ^ m=-“" n=-«> V. ^ 






+ 2(m - n), - — + 2(m + n) 


-2].Kw|- 


tan X, - K + 2(m - n), - + 2(m + n) 


■3 

Ky^, “ ^ + 2(m + n) - 1,^ - 2(ni - n) - ij + “ tan - ^ + 2(m + n) - 1,^ - 2(m - n^ 

Kyj^,^ - 2(m + n) + 1, - ^ + 2{m - n) - ^ - Ky^ - tan x,^ - 2(m + n) + 1, - ^ + 2(m - n)j 


-2tKu 


lx = 9QO’ 


tan X>^ - 2(m ~ ~ 2(m + n) + l| + 2tKy 


■] 


: 90® 


^ - tan x,g - 2(m + n) + 1. 


- i + 2(m - n)j 


(7) 
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where for 3-^® from equations (1) and (3) and t = 0. For 

and Ky are from equations (4) and (6) and t = 1. 


/ OO OO 

- ").§ - * "3 

\m=0 n=0 



+ n) + ^ 


S^m=n?^0 


Kuji, - i + 2(m - n), - ^ + 2(m + n) - ^ - Ku^ - tan x, - ^ + 2(m - n), - ^ + 2(m + n) - ij 
• Ku^, - ^ + 2(m + n) - 1,^ - 2(m - n) - ij - Ku^ - tan x, - ^ + 2{m + n) - 1,^ - 2(m - n)j 
■ Ku^,^ - 2(m + n) + 1, - ^ + 2(m - n) - ij - Kuj|j - tan x,^ - 2(m + n) + 1, - ^ + 2(m - n^ 


■ 4tKu 


lx = 90° 


^ - tan x,^ - 2(m - n),^ - 2(m + n) + iju 


( 8 ) 


where, for 5u l, K^ is from equations (2) and t = 0, and, where, for 6u,D> ^u is 
from equations (5) and t = 1. 

When the wake passes directly rearward (x = 90®), all the terms of equations (7) 
and (8) in which the x-dimension appears as ^ - tan x are zero and may be ignored in 

numerical calculations. 

When m = n = 0, the first term represents the wake in free air. Observe that this 
term is omitted in equations (7) and (8) since it is desired that these equations represent 
only the additional interference caused by the walls. 

Equations (7) and (8) represent the interference only in a closed tunnel. Equivalent 
expressions could be obtained for an open tunnel; however, as discussed in references 1 
and 3, the assumptions inherent in the theory for that case would make the results of 
dubious validity. 


Circular Tunnel 

To approximate the circular tunnel, the interferences are assumed to be equal to the 
average interferences in square and diamond tunnels of equal cross-sectional area. 

For equal areas, ttR^ = 4h 2 = 2D^ so that 

H 

R 2 i 

> (9) 



Thus, for any given value of x/R, y/R, and z/R in the circular tunnel, equa- 
tions (9) may be used to convert the given coordinates for proper entry into the equations 
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of reference 3 or 5 for the square tunnel or into equations (7) and (8) of the present paper 
for the diamond tunnel. After calculating both sets of interference factors, the average 
values are assumed to represent the interference factors for the circular tunnel. A 
FORTRAN program for performing the calculations is presented in the appendix. 

RESULTS AND DISCUSSION 
Interference at the Model 

Figure 3 compares the interference factors at the model for all three tunnels. At 
the higher wake angles, the interference factors for vertical interference are somewhat 
larger in magnitude in the round tunnel than in the square tunnel; whereas, the reverse 
is true at very low skew angles. This result is in accordance with the previously antic- 
ipated trends. The factor for horizontal interference due to lift 6^^L always some- 
what less in the round tunnel and the factor for horizontal interference due to drag 6 u,d 
is always more negative in the round tunnel than in the square tunnel. 

Two considerations should be noted carefully when the interference factors pre- 
sented in figure 3 are examined. First, an effective wake angle should be used in applying 
the results in order to account for the effect of wake rollup on the inclination of the wake 
(refs. 13 and 14), and this effective wake angle results in a substantial increase in y. 
Secondly, it has been shown experimentally (refs. 11 and 15) and theoretically (ref. 12) 
that there is a minimum value of y, which can be tolerated in V/STOL testing, beyond 
which the flow in the tunnel is so distorted that the measured data become meaningless. 
Thus, in a practical sense, it is only necessary to consider the reasonably large values 
of y. Under these circumstances, examination of figure 3 indicates that the interference 
in the round tunnel will not be greatly different from that in the square tunnel. For most 
practical cases, the interference factors will differ by only about 10 percent between the 
two tunnels. Thus, if the model is reasonably small, with a span on the order of one- 
fourth or one-third of the tunnel diameter, it should be completely satisfactory to use the 
interference factors for a square tunnel. Indeed, this procedure is quite attractive since 
detailed computer programs applicable to the square tunnel are available. (See refs. 4 
and 5.) 

At a wake angle of 90°, the present result for the square tunnel is 5-w,L = -0.545 
and for the diamond tunnel is 6-^^, L = -0.547. Because of the different definitions of the 
interference factors, it is necessary to divide these values by -4 (ref. 3) before comparing 
them with Glauert's value of 0.137 for the square tunnel. The aforementioned division 
leads to satisfactory agreement with Glauert's result when the differing degrees of 
numerical approximation are considered. It is interesting to note that there is no sig- 
nificant difference between the interferences in the square and diamond sections when the 
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wake is undeflected (x = 90°). In either case, the interference obtained is within about 
10 percent of Prandtl's (ref. 6) exact value of 0.125 for the circular tunnel. 

Interference Along Longitudinal Axis 

Similar comparisons are shown for the longitudinal distribution of interference 
along the tunnel axis in figures 4 to 7. Not all the computed points are shown; the curves 
are drawn from calculations made at 29 values of x/R. For typical tail lengths, on the 
order of 1/2 R to R, examination of these figures indicates that corrections at the tail 
for pitching moment may be somewhat less satisfactory than those at the center of lift if 
the factors for a square tunnel are applied to tests in a circular tunnel. The error will 
be reduced, however, as the wake angle increases toward 90*^, and such correction is 
obviously superior to no correction at all. 

Interference Along Lateral Axis 

Eqviivalent comparisons on the lateral axis of the tunnel are shown in figures 8 to 11. 
The curves are drawn from calculations made at 13 values of y/R. In general, the trends 
of the distribution over this axis of the factors related to longitudinal interference are 
about the same; the major difference being in the overall level of interference. Unfortu- 
nately, the interference in the square and diamond test sections tends to produce opposite 
trends in the lateral distribution of the interference factors related to the vertical inter- 
ference velocity; the differences become greater as the wake skew angle increases. The 
approximate solution for the circular tunnel is rather uniform because it is taken merely 
as the average of the square and diamond sections. The sensitivity of the result to the 
sidewall configuration leaves some doubt that similar uniform results would be obtained 
for a finite-span system, the representation of which would require laterally offset ele- 
mental wakes (ref. 4). Thus, the present results should be construed as applying only to 
lifting systems of small lateral extent. 

CONCLUDING REMARKS 

An approximate treatment of V/STOL wall interference in a circular tunnel indi- 
cates that the interference factors at the model for this tunnel should be of the same 
order of magnitude as those presently available for a square tunnel of equal cross- 
sectional area. There is a greater degree of uncertainty with respect to the lateral and 
longitudinal distributions of interference; however, the available results for the square 
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tunnel should be reasonably close to those of the circular tunnel provided that the model 
is relatively small in comparison with the test section. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 31, 1970. 
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APPENDIX 


FORTRAN PROGRAM FOR APPROXIMATING THE WIND-TUNNEL INTERFERENCE 
FOR A SMALL MODEL IN A CIRCULAR TUNNEL 


THIS PRCGkA^ hAS WRITTEN IN CCC FQRTRANt VERSION 2. It TO RUN CN CDC 60C0 
SERIES COMPUTERS WITH TEE SCOPE 3.0 OPERATING SYSTEM AND LIBRARY TAPE. MINOR 
EODIEICATIUNS MAY BE PECUIREO PRIOR TO USE IN OTHER COMPUTERS. THIS PROGRAM 
HAS BEEN EUUNC SATISFACTORY ON THE AFOREMENTIONED COMPUTERS WHICE CARRY THE 
ECUIVALENT OF A F F PC X IM A T ELY 15 DECIMAL DIGITS. COMPUTERS OF LESSER PRECISION 
MAY kEUUIRE MCDIFICATICN TO DOUBLE PRECISION IN ORDER TO OBTAIN RESULTS OF EQUAL 
NUMERICAL ACCURACY. 

TWO SUBROUTINES ARE USED. SUBROUTINE <SCUARE> IS A MODIFICATION OF THE 
FROGRAM GIVEN AS APPENDIXES A AND C OF NASA TECHNICAL MEMURANOUM X-1740. IT 
CCMPUTES THE INTERFERENCE FACTORS FCR A SQUARE, CLOSED TEST SECTION. SUBROUTINE 
<CIAMCNU> IS A FORTRAN COCING OF EQUATIONS (71 AND (8) OF THE PRESENT PAPER. 
CONVERSIONS BETWEEN NCN C I ME NS ION A L I / A T I ON S WITH RESPECT TO R, H, AND D ARE 
MADE IN THE APPROPRIATE ROUTINES. 

THE REQUIRED INPUT, FOUNC AT ADDRESS 1 (LINE (A 9)) IN FORMAT 103 (LINE 
(A CONSISTS CNLY OF THE X- , Y-, AND Z-DISTANCES FROM THE MODEL TO THE 

FLINT AT rtHICF TFE INTERFERENCE FACTORS ARE REQUIRED. ALL THREE DISTANCES 
MLST BE NCND IMLNS IGNAl. I ZEC WITH RESPECT TO THE RADIUS OF THE CIRCULAR TEST 
SECTION. 


PROGRAM RLLNU 1 I N P L T ,U LTP U T , T AP E5= I NPU T , T AP E6 = 0UTPUT ) (A 1) 

C (A 2) 

C (A3) 

CLMMCN XLVERH,yOVERH,ZOVERF,DELTA(A) (A 

UIMENSICN CEL I A1 (A ) ,DELTA2 ( A ),NAMEI A) , C ( 11 ) (A 5) 

DATA ( NAM £ ( 1 ) , 1=1 ,A )/lCFDFLTA ( W ,L) , lOHDELTAIU ,L ) , lOhDELTAI W,C ) , lOF (A 6) 

IDELTAIU ,L )/ (A 7) 

DATA ( C ( 1 ) , 1= 1 , 11 )/-3. ,3. , 10 .,20. , 30. , AC. , 50. ,6C. , 70. ,80., 9C. / (A 8) 

I READ (5,1C3) XCVERR, YOVERR ,Z0VERR (A 9) 

IF (EOF, 5) 99S,2 (A 10) 

API=3.1^155Zt£2!B979 (All) 

HCVF.RR = SLRT ( FI/A.) (A 12) 

XCV£RH=XUV£RR/FCVERR (A 13) 

YUVERh=YLV£RR/FCVERR (A lA) 

ZUV£RH=ZO VERR/FCVERF (A 15) 

WRITE (t,M8) XOVERR , YLVERR,ZaVERR (A 16) 

WRITE (6, 15C ) ( A 17) 

DO A1 K=l,ll (A 18) 

WRITE (6,1A9) C(K) (A 19) 

DO 20 1=1,9 (A 20) 

20 DEL TA( I )=DELTA1( 1 ) = DELTA2( I ) =0. (A 21) 

C (A 22) 

C ( A 23) 

CALL SCUAPE (C(K ) ) (A 2A) 

C ( A 25) 

C (A 26) 

DU 3 1=1 ,A (A 27) 

3 DELTA1(I)=CELTA(I ) (A 28) 

C (A 29) 

C (A 30) 

CALL DIAMCNC (CIK)) (A 31) 

C (A 32) 
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APPENDIX — Continued 


C (A ?3) 

UC ‘i 1=1,4 (A 34) 

UELTA2( I»=CfcLTA(I ) (A 35) 

4 UELTA( I )= (CELTAK I )*OELTA( I ) )/2. (A 36) 

WHITE (6,151) (( NAHE( n ,DELTA1 (I ) ,UELTA2( I) ,DELTA( I ) ) , 1=1 ,4 ) (A 37) 

41 CUNTINtE (A 38) 

GC TC 1 (A 39) 

1U3 t-OKHAT (3F1C.3 ) (A 40) 

148 FCH^AT (1M//22X*INTERFEKENCE FACTUKS AT A POINT NEAR A VAMSHINGL (A 41) 

lY SPALL PCC£L<'//32X*APPK0XIPATE SCLOTICN IN A KCUNO CLOSED TUNNEL* (A 42) 
2 ///23X*X/K =*F8.3,9X*Y/R =*F 8 . 3 , 1 1 X* Z /R =*F8,3//) (A 43) 

149 FCRPAT (//EX+CFI =*F7.3/) (A 44) 

130 FCRPAT (22X*SCLAKE TONNE L ♦ 1 9 X*0 I A MC NO T UNN EL* 1 5X* AV ER A GE (RCLND TU (A 45) 

2NNEL)*/2<cX,13 (1F=),19X,14( IF=),13X,22< 1H=)) (A 46) 

131 FORMAT (10,A1C,F12.4,2F32.4) (A 47) 

‘>99 STOP (A 48) 

tNO (A 49) 


SUdKOUTINE CIAPCNC (ANGL) (A 50) 

CCPPCN XCVEPF,YCVERF,Z0VERF,DELTA(4) (A 51) 

DIMENSION V(4) (A 52) 

SC=SIN (ANGL*0. 0174532925199 ) (A 53) 

CC=COS ( ANGl*0.C174532923199 ) (A 54) 

XUVERO=XCiVERF/£QR7(2.0) (A 55) 

Y0VER0=Y(jVERF/SCRT(2.0 ) (A 56) 

ZOVERC=ZCVEPF/SCRT(2.0) (A 57) 

OU 1 1=1,28 (A 58) 

1 DELTA( I )=C .C (A 59) 

00 2 IM=1,7 (A 60) 

DO 2IN=1,7 (A61) 

DO 3 1=1,4 (A 62) 

3 V(I )=0.0 (A 63) 

M = 2.*( IP- IN ) (A 64) 

N=2.*( IN+IM-8 ) (A 65) 

IF (IM.EU.4.ANC.IN.EQ.4) GO TO 4 (A 66) 

X=XCVERC (A 67) 

V = YOV£RL--P (A 68) 

Z=ZLVERL'-N (A 69) 

A = SCKT (X*X-#Y*Y*Z*Z ) (A 70) 

B=A+Z*CC-X*£C (A 71) 

V ( 1 )= ( ( X*>4 t*Y ) /( B*A*A*A) )- { (Z+A*CC )/ ( B*A) ) **2 (A 72) 

V (2 )=V (3) =- (X*Z )/ ( e*A*A*A )- (Z+A*CC ) ♦( X-A*SC )/ ( B*B*A*A) (A 73) 

V(4)=( ( Y*y ^Z*Z )/( B*A*A*A) )-( (X-A*SC)/ ( e*A) ) **2 (A 74) 

DC 5 1=1,4 (A 75) 

3 DELTA! I )=CELTA ( I )+V( I ) (A 76) 

4 IF ( ANGL.EC .9C . ) GC TO 6 (A 77) 

X = XGVERC- (SC/CC ) (A 78) 

Y=YCVERD-P (A 79) 

Z = ZCVEPC-N41 . (A 80) 

A=SQRT(X*X+Y*Y+Z*Z) (A 81) 

B=A+Z*CG-X *SC (A 82) 

V (1 ) = ( ( X*X + y*Y )/( e*A*A*A) )- ( (Z + A*CC)/ ( B*A ) )**2 (A 83) 

V (2 ) = V ( 3)=- ( X*2 )/ (8*A*A*A)-( Z+A*CC ) *( X-A*SC) / ( B*B*A*A) (A 84) 

V (4 ) = ( ( Y*Y4Z*Z )/( e*A*A*A) )- ( (X-A*SC )/ ( B*A ) )**2 (A 85) 

CO 7 1=1,4 (A 86) 

7 DELTA! I ) = C£LTA(I)-V(I) !A 87) 

fa X=XCVERC !A 88) 

Y=-YOVEPU+P !A 89) 

Z=-ZCVERD*N-2 . !A 90) 

A=SOKT ! X* > <y*Y4Z*Z ) !A 91) 

6=A+Z*CC-X*SC !A 92) 
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APPENDIX — Continued 


V(l) = ( (>*>4>f*Y)/(B*/5*A*A))-((Z + A*CC)/(e»A))»*2 
V(2)=V(3) = -(X4Z )/(E+A*A«A)-(Z+A*CC )♦( X-A^ SC) / ( B*B*A* A ) 

V (4)=( (V*V4Z^Z )/( B^A*A«A) )-( (X-A+SC)/ ( B^A) )**2 
DO 8 i=l,^ 

6 DELTA ( I )=CELTA( !)♦« (-!• )=»♦! )*V( I) 

IF (ANGL.EC.9C.) GC TO 9 
X=XCVERC- (SC/CC ) 

Z=Z4l. 

A=SURT(X*X4Y*V4Z*Z) 

B=A+Z*CC-X’»SC 

V(U = ( (X*>4V*V)/( E*A+A*A) )-( (Z4A*CC)/(B«A) )**2 

V (2 ) = V ( 3) = - (X<Z ) / (e*A*A«A)-(Z+A«CC )♦( X-A’i'SC ) / ( B*B*A^A ) 

V (4)=( {y*Y4Z*Z )/{ B^A^A*A) )-( (X-A*SC)/ ( B*A) )**2 
DC 10 I=it^ 

10 DELTA ( 1 )=CtLT/> ( I)-< (~1 . n 

9 X=XCVERC 

Y=-Z0VERC4^-^ . 

Z=YGV£PC-)'-l. 

A^SGKT (X^X4Y*Y4Z^Z) 

8=A+Z*CC-X*SC 

V( i )=-Y*( Z4A+CC )/ (A#A+B^8)-Y^Z/(B*A*A*A) 

V(2 )=-(X*Z )/(6#A«A*A)-( Z+A*CC)*(X-A*SC )/(B*B*A*A| 

V(3 )=-Y*( X-A+SC )/ (AtA*B*8)-Y^X/(B*A*A*A) 

V(4 )^( (Y+Y4Z+Z )/( B’^A^A + A) )-( ( X- A*SC ) / ( 8 ♦A ) )^*2 
00 11 1=1,4 

11 DELTA! 1 )=CELT/! ( n + ( (-1. )**I)*V( I) 

IF (ANGL- EC.9C.) GC TO 12 
X=XCVERL- (SC/CC ) 

Z=Z+1* 

A=SGKT ( X*X4Y*Y4Z*Z ) 

B=A+Z*CC-X*SC 

V( I )=-Y*( Z+A*CC )/ ( A*A^B*B)-Y*Z/ (8*A*A^A) 

V(2) = -(X<^Z )/(E*A4'A*A)-(Z+A*CC)*(X-A*SC)/(B*8*A*A) 

V(3) = -Y*( X-A«SC )/ ( A*A*B*B)-Y*X/ (B+A*A*A ) 

V(4) = ( (Y*Y4Z + Z)/(B*A*A*A))-( { X- A=^ SC) / ( B *A ) ) ♦♦2 
DO li 1=1,4 

13 DELTA! 1 ) = C E 1 T A ( I) - ( (- 1 . J «♦ n ( I) 

12 X=XQVERC 
Y=ZCVERC-A41. 

Z=-YUVEPD4F-1 . 

A=SCRT (X*X 4Y*Y4Z*Z ) 

B=A+Z^CC-X*SC 

V! 1 )=-Y*( Z4A*CC )/( A^A+B*B)-Y*Z/!B*A*A«A) 

V (2 )=- (X^Z )/ ( E<A#A^A)-(Z4A*CC)*(X-A*SC )/! B*B*A«A) 

V(3 )=-Y*( X-A«SC )/(A^A*B«0)-Y*X/(B4^A*A*A) 

V(4)=( (Y*Y4Z*Z)/!e^A«A*A))-( (X-A+SC)/! B+A) )**2 
DO 14 1=1,4 

14 DELTA! I )=CELTA! I)4V( I ) 

IF (ANGL.EC.9C.) GC TO 2 
X=XOVERC-(SC/CC) 

Z=Z+1. 

A=SCRT !X*>4Y*Y4Z*Z) 

B=A+Z*CC“X*5C 

V(l )=-Y^( Z4A*CC)/(A*A^B*0)-Y*Z/ (B^A«A+A) 

V!2)=-(X<=Z )/(E«A^A*A)-( Z4A*CC)*(X-A*SC)/(B«6*A*A) 

V(3 )=-Y*!X-A^SC )/ (A*A*8*B)-Y#X/(B*A^A*A) 

V(4) = ( !Y<^Y4Z»Z)/!B*A*A*A))-( (X-A+SC )/ ( B*A ) )*«2 
DC lb 1=1,4 

lb DELTA! 1 )=CEL1/ ( n-V( I ) 

Y=Y0VERC-F 

Z=ZLVERC-N41, 

A=SCKT !X«X4Y*Y4Z+Z) 


(A 93) 
!A 94) 
(A 95) 
!A 96) 
(A 97) 
!A 98) 
(A 99) 
(A 100) 
!A 101) 
(A 102) 
(A 103) 
(A 104) 
(A 105) 
(A 106) 
(A 107) 
(A 108) 
!A 109) 
!A 110) 
!A 111) 
(A 112) 
(A 113) 
!A 114) 
(A 115) 
(A 116) 
(A 117) 
!A 118) 
(A 119) 
(A 120) 
(A 121) 
(A 122) 
(A 123) 
(A 124) 
(A 125) 
(A 126) 
(A 127) 
(A 128) 
!A 129) 
(A 130) 
(A 131) 
(A 132) 
(A 133) 
(A 134) 
!A 135) 
(A 136) 
(A 137) 
(A 138) 
<A 139) 
(A 140) 
(A 141) 
(A 142) 
(A 143) 
(A 144) 
(A 145) 
(A 146) 
!A 147) 
(A 148) 
(A 149) 
(A 150) 
(A 151) 
(A 152) 
(A 153) 
(A 154) 
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APPENDIX — Continued 


B=A-X (A 155) 

V (3 )=Z/(A*A^A ) (A 156) 

V(4)=X/(A«A*A ) (A 157) 

CC 16 I=3»4 CA 158) 

lo DELTA( n=CtlTA( r)+4.^V( I) <A 159) 

2 CGNTlNbt (A 160) 

UG 19 1=1,4 (A 161) 

19 UELTA( I )=-CELTA(l )/3. 14 159265358979 (A 162) 

RETURN $ENC CA 163) 


SUBRUUT INE SQUARE (ANGL) 

CCMMON XU V ERE , YC V ERE , ZO VERH , CELT A ( 4 ) 

DIMENSICN V(3,9) ,ACEL( 4) 

SC=SIN( ANGL*0*C174532925199) 

CC=COS( ANGL *0.0174532925199) 

Z6= ZCVERE41. 

Ze=-Z6 
Z7 = Z8-1. 

DC 8 Jl=l,4 
8 DELTA! J1)=C. 

DG 10 M=l,7 
DO 10 N = 1 ,7 

IF (N.EC.4.ANC.M.EQ.4) GO TO 10 
DC 11 Jl=l,3 
DO 11 J2=l,9 

11 V(J1,J2)=0. 

DO 12 Jl=l,4 

12 ADELCJ 1)=C. 

AK=H-4 

AN=N-4 

X= XOVERH 

Y= YCVERH-2.*AE 

Z= (ZGVERH-4.*AN) 

A=SCRT (X*X4Y*Y*Z*2 ) 

B=A+Z*CC-X*SC 

VU f 1 ) = ( C X*X4Y*Y )/( B*A*A*A) )-( C Z+A*CC) /(8*A ) ) **2 
V(2,l) =-< X*Z )/(B*A*A*A)-(Z + A*CC)*(X-A*SC)/(B*B*A*A). 
V(3,l )=(( Y*Y+Z*Z)/( e*A*A*A) )-( ( X-A*SC )/(B*A) )**2 
Z=-Z-2 . 

A=SQRT (X*X4Y*Y4Z*Z ) 

8=A+Z*CC-X*SC 

V(l,3) = ((X*X + Y*Y)/Ce*A*A*A) )- C { Z4A*CC ) / ( B*A ) ) **2 
V(2,3)=-( X*Z )/(e*A*A*A)-(Z+A*CC)*(X-A*SC)/(B*B*A*A) 
V(3,3)=((Y*Y4Z*Z)/(E*A*A*A))-((X-A*SC)/(B*A))**2 
IF (ANGL.EC.9C.C) GO TO 13 
X=X-(SC/CC ) 

z=-z-i. 

A=SCRT(X*X4Y*Y+Z*Z) 

B=A+Z*CC-X*SC 

V(l,2)=((>*>4Y*Y)/(e*A*A*A))“(( Z+A*CC )/ (b*A ) )**2 
V(2,2)=-(X*Z)/lB*A*A*A)-(Z+A*CC)*(X-A*SC)/{ B*B*A*A) 
V(3,2) = (( Y*Y4Z*Z )/( e*A*A*A) )-( ( X- A * SC ) / ( t3* A ) )**2 
B=A-X 

V(l,5) = (( X*X4Y*Y)/(0*A*A*A) )-(Z/(B*A) ) **2 
V(2,5)=Z/( A*A<A) 

V(3,5)=X/ ( A*A*A) 

Z=-Z 

B= A+Z*CC“X*SC 

VU,4) = ((X*> + Y*Y)/(e*A*A*A))-“((Z + A*CC)/(8*A))**2 
V(2,4)=-( X*Z )/(B*A*A*A)-'(Z + A*CC)*(X-A*SC)/(B*8*A*A) 


(A 164) 
(A 165) 
(A 166) 
CA 167) 
CA 168) 
CA 169) 
CA 170) 
CA 171) 
CA 172) 
CA 173) 
CA 174) 
CA 175) 
CA 176) 
CA 177) 
(A 178) 
CA 179) 
(A 180) 
CA 181) 
CA 182) 
CA 183) 
(A 184) 
CA 185) 
CA 186) 
(A 187) 
CA 188) 
CA 189) 
CA ISO) 
(A 191) 
(A 192) 
CA 193) 
CA 194) 
CA 195) 
CA 196) 
CA 197) 
CA 198) 
CA 199) 
(A 2C0) 
CA 201) 
(A 202) 
CA 203) 
(A 2C4) 
CA 205) 
CA 206) 
(A 207) 
CA 208) 
(A 209) 
CA 210) 
CA 211) 
(A 212) 
(A 213) 
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APPENDIX — Concluded 


V<3,4)=((Y*Y+2*Z)/(e^A*A«A) )-( (X-A^SC )/(B«A) )**2 

13 ADEL(1)^V( 1,U-VC 1 ,2 V ( I , 3 ) ♦ V t 1 , 4 ) 

ADEL(2 )=V(2f l)-VC 2,2)+V(2» 3)-V( 2,4) 

ADEL(3)=V (2,1 )-V( 2,2)-V (2,3) +V( 2f 4) *2.*V(2, 5) 
ACEL(4) = VO,l )-V(3,2)+V(3,3)-Vl 3,4)+2.«V(3,5> 

00 14 Jl=l,4 

14 DELTAC Jl)=CELTA(Jl)+ADEL( Jl) 

10 CONTINUE 

X= XOYERH 

Y= YGVERH 

Z=Z7 

A=^SCRT(X*X4Y4Y4Z*ZI 

B=A+Z^CC-X^SC 

VI 1,7J = (( M-X + Y^Y )/(B*A*A*A))-( ( Z+A*CC) / ( 8*A )) **2 
V(2,7)=-(X42 ) / (e*A*A^A)-(Z + A«CC )♦ ( X-A^SC ) / ( B*B*A<^A ) 
V(3,7) = ( C Y*Y4 2^Z)/( e^A+A*A) )- ( I X-A*SC ) / ( 8*A ) )**2 
IF (ANGL*EC.9C«C) GO TO 16 
X^X-(SC/CC ) 

Z = Z6 

A=SGRT (X^X4Y*Y + Z + Z ) 

B=A+Z^CC~X«SC 

V ( 1,6)= ( ( XO + Y4Y )/<e4«A*A*A) )-( ( Z+A^CC )/ (B*A) )**2 
VC2,6)=-CX*Z )/(e«A^A*A)-(Z+A4CC)*(X-A*SC)/(8«B<'A*A) 
V(3,6) = ((Y*>4Z*Z )/(e*A^A*A) ) - ( ( X- A« SC ) / { B^ A ) )* + 2 
B=A~X 

V(l,9) = ((>*> + Y*Y)/(etA«A*A))--(Z/(B*A))*42 
V(2,9)=Z/ (/^A*A) 

V(3,9)=X/ (Z + A^/) 

Z=Z8 

B=A+Z*CC~X^SC 

V(1 ,8) =(( >*> + Y^Y I/ie^A’frA + A) )-(( Z + A*CC)/(B*A) )*«2 
V(2,8) = -|X«Z )/(B*A*A*A)-IZ+A*CC X~A*SC ) / ( B«B* A* A) 

V(3,8)=((Y^Y4Z^Z )/(B«A*A4«A) )-((X-A«SC)/(8*A) )**2 

16 ADEL(l)=*-V(l,t)-V(i,7)4V(l,8) 

ADEL (2 )=-V(2,6 )4V(2,7)-V(2, 8) 

ADEL( 3)=^V (2,6)-V(2,7) 4V(2 ,8)42,*V(2,9) 

ADEL (4 )=-V (3,6 )4V( 3,7 )~V( 3,8)42.*V( 3,9) 

00 17 Jl=l,4 

17 DELTA( Jl) = CELTA( Jl )4ADEL( Jl) 

AMT=-2. /3. 14 159266358979 

DO 19 Jl=l,4 

19 DELTA! Jl) = A^<T«CELTA(Jl) 

RETURN 

END 


(A 214) 
(A 215) 
(A 216) 
(A 217) 
(A 218) 
(A 219) 
(A 220) 
(A 221) 
(A 222) 
(A 223) 
(A 224) 
(A 225) 
(A 226) 
(A 227) 
(A 228) 
(A 229) 
(A 230) 
(A 231) 
(A 232) 
(A 233) 
(A 234) 
(A 235) 
(A 236) 
(A 237) 
(A 238) 
(A 239) 
(A 240) 
(A 241) 
(A 242) 
(A 243) 
(A 244) 
(A 245) 
(A 246) 
(A 247) 
(A 248) 
(A 249) 
(A 250) 
(A 251) 
(A 252 ) 
(A 253) 
(A 254) 
(A 255) 
(A 256) 
(A 257) 
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(b) Wake and image in ground effect. 

Figure 2.- Wake and image systems used in developing interference factors 

for a diamond test section. 
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(c) Vertical interference due to drag. 
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(d) Horizontal interference due to drag. 


Figure 3.- Interference factors at model centered in tunnel. 
Symbols show the shape of test section. 
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Figure 8.- Distribution of interference factors over lateral axis of tunnel. Shape of symbols denotes 
shape of test section, x = (Note that for this skew angle, 6 u,l = 0-) 
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Figure 11.- Distribution of interference factors along lateral axis of tunnel. Shape of 
symbols denotes shape of test section, x = x = 6u,L> ^w,D) 

6u,D 3^re all zero over the entire axis.) 
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